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Abstract 

We present a measurement of the branching fraction and time-dependent CP violation param- 
eters for B'^ — > D*~^ D*~ Kg decays. These results are obtained from a 414 fb~^ data sample that 
contains 449 x 10^ BB pairs collected at the T (45*) resonance with the Belle detector at the KEKB 
asymmetric-energy e'^e~ collider. We obtain the branching fraction, 

13 (B^ D*+D*-K^s) = [3-4 ± 0.4 (stat) ± 0.7 (syst)] x lO'^, 

which is in agreement with the current world average. We also obtain an upper limit on the product 
branching fraction for a possible two-body decay, 

B{B^ D+ {2536)D*-)B{D+ (2536) ^ D*+ K^) < 7.1 x 10"^ (90% CL). 

In the traditional 2-parameter time-dependent CP analysis, we measure the CP violation param- 
eters, 

Acp = -O.Oltoil (stat) ± 0.09 (syst) 
Z?sin2(/)i = 0.06l[!;^^ (stat) ± 0.06 (syst). 

No evidence for either mixing-induced or direct CP violation is found. In a 3-parameter fit sensitive 
to cos 2^1 performed in the half-Dalitz spaces, s~ < and s~ > where = {D*^Kg), 
we extract the CP violation parameters, 

Jc/Jo = O.mtlfs (stat) ± 0.08 (syst) 
2Jsi/Josin2(/>i = -0.17l[!:^^ (stat) ±0.09 (syst) 
2Js2/Jo cos 2(^1 = -0.23;[|;^f (stat) ±0.13 (syst). 

A large value of Jc/Jo would indicate a significant resonant contribution from a broad unknown 
D**~^ state. Although the sign of the factor, 2Js2/Jo, can be deduced from theory, no conclusion 
can be drawn regarding the sign of cos 2(f)i given the errors. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 
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CP violation in the Standard Model is due to a complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix [l|, At present, CP violation has been clearly 
observed by the BaBar [3,] and Belle j3] collaborations in the decay — > J/ip K^, while 
many other modes provide additional information on CP violating parameters. One such 
decay is D*^ D*~ K^, where sin20i can, in principle, be extracted from the time- 

dependent rate asymmetry |5|, 

r ft) - r (t) 

p ^ ^ = Acp COS (Amdt) + D sin 20i sin (Amat). (1) 

r (t) is the decay rate for B^ D*^ D*^ Kg at proper time, t, after production while T {t) 
represents the charge conjugate rate. Amd is the mass difference between the two B^ mass 
eigenstates, Acp is the direct CP violating component and sin 201 measures mixing induced 
CP violation. The factor, D, is the dilution of the CP asymmetry and may arise from two 
sources. The contribution from different waves in the amplitude of B^ D*~^D*^Kg 
influences the CP mixture and polarization of the final state. Furthermore, if an inter- 
mediate D*+Kl, resonance exists, additional dilution of the CP asymmetry occurs as the 
decay becomes self-tagging. CP contamination from penguins and final state interactions 
are expected to be small in this mode jl]. 

If resonant structure were to exist, cos20i may also be extracted, assuming no direct CP 
violation, from the time-dependent rate asymmetry [5|, 



T(t)-T(t) Jc , /2Jsi , 2J,2 , \ X 

f,,,. , J,,.. =Vy-r cos {Arridt) - —— sm20i + Vy^— cos20i sm (Amat), (2) 

i [t) + L [t) Jq V Jq -'o / 

in the half-Dalitz space, s~ < s"*" or > s+. The Dalitz variables are defined as = 
{D*'^Kg) and rjy = +1 (—1) if < (s^ > s"^). Jc, Jq, Jsi and Js2 are integrals over 
the half-Dalitz space, > s^, of |ap + |ap, |ap — |ap, the real component, 3fJ(aa*) and 
the imaginary component, Q{aa*), respectively, where a (a) are the decay amplitudes of 
B° (5°) D*+D*~K^s- Note that J^2 is non-zero only if 5° D*+D*~K^s has a resonant 
component: in this case, Jc may be large. 

This measurement of the branching fraction and CP violating parameters in B^ — »• 
D*+D*-X^ is based on a 414 fb"^ data sample that contains 449 x 10^ BB pairs, collected 
with the Belle detector at the KEKB asymmetric-energy e+e" (3.5 on 8 GeV) coUider Q. 
Operating with a peak luminosity that exceeds 1.6 x 10^^ cm~^s~^, the collider produces the 
T (45) resonance (a/s = 10.58 GeV) with a Lorentz boost of /?7 = 0.425, opposite to the 
positron beamline direction, z. Since the B^ and B^ mesons are approximately at rest in the 
T (45) Center-of-Mass System (CMS), the difference in decay time between the BB pair. 
At, can be determined from the displacement in z between the final state decay vertices. 

At ^ ("^^ ~ "^-^^ ^ ^. (3) 

The subscripts, CP and Tag, denote the reconstructed final state, Bcp, and the final state 
from which the flavor is identified, -Bxag- The quantity. At, can be substituted for t in ([T]) 
and ([2]) due to coherent BB production at the T {AS) resonance [6|. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
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(TOF), and an electromagnetic calorimeter comprised of Csl (Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [sl]. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector (SVDl) was used for 
the first sample of 152 x 10^55 pairs, while a 1.5 cm beampipe, a 4-layer silicon detector 
(SVD2) and a small-cell inner drift chamber was used to record the remaining 297 xlO^-B-B 
pairs ^]. 

We reconstruct D*^D*^Kg from D*+ D^n^^^ and D'^n^i^^ requiring at 

least one D^n^^^ decay. and candidates are reconstructed from K~tc^, 
/r^TT+TT", i^~7r+7r°, K^TT^TT^TT^ , K" and K~7r~^7r~^, K~K^ vr^; cases with two 

— >■ Kgir^n' decays are rejected. We also reconstruct Kg —>■ tt+tt^ and tt" 77. The 
reconstruction of the charge conjugates in the above decay chain is implied. 

According to Monte Carlo (MC) studies, continuum background does not give a significant 
contribution to this mode, so only a loose requirement on the ratio of the second to zeroth 



Fox- Wolfram moment [10| is used, R2 < 0.4. Charged tracks satisfy loose criteria on their 
impact parameters relative to the interaction point (IP), \dr\ < 0.4 cm and \dz\ < 5.0 cm 
with some additional SVD requirements [ll|. With information obtained from the CDC, 
ACC and TOF, particle identification (PID), or K/ir separation is determined with the 
likelihood ratio, Ck/{Ck + C,,). Here, £j is the likelihood that the particle is of type i. 
Candidate vr'^'s are selected from pairs of photons with energies greater than 30 MeV. The 
7r° momentum is required to be greater than 200 MeV/c in the laboratory frame. Charged 
slow pions used in the reconstruction of i5*'s, are not subject to impact parameter or PID 
cuts. Similarly, the vr" momentum cut is not applied to slow 7r°'s used to reconstruct _D*'s. 
We generate signal MC with EvtGen l2] using a detector simulation based on 



GEANT3 [13|. All mass windows are chosen from correctly reconstructed signal MC by 
fitting the mass distributions with double Gaussians and applying a 3cr cut on the wider 
Gaussian. This method is applied to tt*^, K^, D and the D* — D mass difference for D* 
selection. When multiple B mesons are reconstructed, the candidate with the smallest value 

XI^^^XHX), wherex-(A-)^ l— (4) 

is selected. We sum over the prompt K^, the two daughter D mesons and the two daughter 
D*'s where a {X) is the width of the narrow Gaussian found in the fit to the invariant mass 
of particle X. 

As we reconstruct the decay chain, vertexing algorithms are applied to and D mesons 
from their charged daughters. The vertex of the B is determined from the pseudo-tracks 
of the two daughter D mesons with an IP constraint. We do not include the slow pions in 
the vertex fit due to their poor resolution. The prompt pseudo-track is also excluded 
as its decay vertex has a relatively large displacement from the IP. However, both the slow 
pion and candidate are constrained to originate from the B vertex determined by the D* 
mesons and their momentum is recalculated to improve the resolution of the B candidate. 

The B decay is described by the variables Mbc = \/(^beam)^ - {Pb^^? and = E"^^^ - 
-^beam- The Mbc sigual shape is modeled by a Gaussian while the background is modeled 



with an ARGUS function In Ai?, the signal shape is represented by a double Gaussian 
and the background with a linear function. With the signal shape determined from MC, 
we define the signal box to be 3cr of the Gaussian in Mbc by 3(7 of the narrow Gaussian in 
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A£^, which is 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and -0.04 GeV < < 0.04 GeV. The 
reconstruction efficiency, which includes all intermediate branching fractions, is determined 
from phase space signal MC to be (6.59 ± 0.13) x 10"^ for SVDl and (9.58 ± 0.16) x 10"^ 
for SVD2, giving an average efficiency of (8.57 ± 0.15) x 10^^. The noticeable difference 
in efficiency between SVDl and SVD2 is due to the better performance of SVD2 with low 
momentum tracking. With the signal shape fixed from signal MC, we perform an extended 
two dimensional maximum likelihood fit to obtain the yield (Fig. [I]), Y = 131.2t\ti (stat) 
events. Using 

e(BO^a.+B.-A-0) = _^, (5) 

where e is the average detection efficiency and N (BB) is the number of BB pairs, we obtain 
the branching fraction, 

B (S° ^ D*+D*-Kl) = [3.4 ± 0.4 (stat) ± 0.7 (syst)] x 10-^ (6) 

which relies on the assumption of equal production of neutral and charged B meson pairs 
from the T (45*) . 




5.24 5.26 5.28 5.3 

M,,(GeV/c^ 




0.04 0.12 0.2 

A E (GeV) 



FIG. 1: Mbc — fit projections. The solid curve shows the fit result. 



Systematic errors are estimated from various sources such as the number of BB events, 



the uncertainty in the efficiency and the PDG errors on branching fractions [15 



The K'^s^ 

ir^, PID and tracking systematics are determined from the efficiency ratio of data to MC as 
calculated by independent studies at Belle. We also include the difference in yield between 
fits with signal shape parameters obtained directly from data, and obtained from MC, as a 
systematic error. The wider A^^ Gaussian in the fit to data is fixed, from MC, relative to 
the narrow floating Gaussian. Peaking background due to fake vr^j^^'s is also accounted for 
by investigating the D*~^ sideband. We also include possible efficiency variations across the 
Dalitz plane by comparing the average efficiency of phase space signal MC with the average 
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efficiency as determined from bins of m {D*'^K^) using data weights. The systematics are 
summarized in Table HI 



Source 


Systematic error (%) 


Number of BB pairs 


1.3 


MC reconstruction efficiency 


1.7 


K^, vr'^and PID efficiency 


5.8 


Daugfiter brandling fractions 


8.4 


Tracking efficiency 


11.9 


Signal shape 


10.4 


vr*^ peaking background 


4.7 


Dalitz efficiency dependence 


2.6 


Total 


19.8 



TABLE I: B^ — > D*^D* Kg branching fraction systematics. 

The D*^Kg spectrum inside the signal box is examined to search for resonant structure. 
Fig. [2](a) shows the background subtracted D*^Kg raw yield versus invariant mass overlaid 
with reconstructed phase space obtained from signal MC. The background is taken from a 
sideband defined to be 5.20 GeV/c^ < Mbc < 5.26 GeV/c^ and -0.2 GeV < AE < 0.2 GeV. 
From this plot, the data inside the signal box does not appear to be entirely consistent with 
phase space, but neither is there any easily identifiable peak. 

We extract the signal yield of Z)^ (2536)D*^ from the mass difference. Am = 

m {D*^ K^)—m {D*^) — m {Kg), in the signal box. Note that while this form of Am has been 
chosen, the Kg was fitted to its nominal mass. A Breit-Wigner function convolved with a 
Gaussian is chosen to model the signal shape, while a simple threshold function, Pskg (Am) = 
cAm^/^, is chosen to describe the background, where c is a normalizing constant. In the fit 
to signal MC, the Breit-Wigner mass and width are fixed to their respective generated values 
while the Gaussian is floated to determine the effects of detector smearing. In the fit to data, 
the same smearing Gaussian is used, but the Breit-Wigner width is conservatively increased 
to the current upper limit on the (2536) width, 2.3 MeV/c^ [l5|. The fitted yield is 
Y (S° (2536)D*-) = 6.2^1^ (stat) events, and is shown in Fig. [SJ^b). The detection 

efficiency in this region is found to be (3.08 ±0.40) x 10"^ for SVDl and (3.98 ±0.36) x 10"^ 
for SVD2, giving an average efficiency of (3.68 ± 0.38) x 10~^, and we obtain the product of 
branching fractions, 

B (5° ^ (2536)D*-)i3 (D+ (2536) ^ D*+ K^g) = [3.8+|l (stat) 1°;? (syst)] x 10-^ (7) 

where the systematics were determined in the same way as for D*^D*^Kg. Thus, 

the upper limit is found to be 

B (fi° ^ D+ {2536)D*-)B {D+ (2536) ^ D*+ i^°) < 7.1 x 10"^ (90% CL). (8) 



To obtain the At distribution, we reconstruct the tag-side vertex from the tracks not 



used to reconstruct the CP-side [ll|, and employ the flavor tagging routine described in 
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FIG. 2: (a) Background subtracted raw yield vs. invariant mass with a superimposed 

reconstructed phase space histogram from signal MC. (b) Measured distribution of m{D*~^Kg) — 
m {D*~^) — m (Kg) in the signal box. The solid line shows the fit results. 



Ref. [la]. Due to imperfect flavor tagging, both a flavor tag, q, and the dilution factor of 
flavor tagging, r, are assigned to each event. The r distribution is divided into 6 bins. 

In the 2-parameter fit, the time- dependent At distribution of D*^ D*^ Kg is given 

by the probability density function (PDF), 

g-|At|/Tgo f 

Vsig (At,g) = <^ l-qAw + q{l-2w)[Acpcos{AmdAt)+Dsm24)ism{AmdAt)] 



Atbo 

(9) 

where q is the tagged flavor, and w, Aw account for imperfect flavor tagging, w is the 
probability of mis-tagging an event. Aw is the difference between and B^ in w and they 
are both determined from flavor specific control samples [l6|. Acp and Dsin20i are the 
two free parameters in this fit. To account for differences between true At and reconstructed 
At, Vsig is convolved with a resolution function described in Ref. 



/-t-oo 
Psig (At', q)Rsig (At - At')diAt'), (10) 
-oo 

which describes the effects of detector resolution and the effect of charm decays on the tag- 
side. The effect of the approximation in Eq. [3]is also included in the resolution function. The 
background PDF is modeled by an exponential lifetime component and a prompt component 
given by a delta function, 

-|At|/TBkg 

PBkg (At) = (1 - fs)— + fs 6 {At - fis). (11) 

^TBkg 

where rekg is the effective lifetime of the background and /i^ sets the At location of the 
prompt component. The PDF, Pskg, is convolved with a resolution function that takes the 
form of a double Gaussian, to give the background PDF, Pskg- The background function 
parameters are determined from the sideband data defined by 5.20 GeV/c^ < Mbc < 5.26 
GeV/c^ and —0.2 GeV < AE < 0.2 GeV. To account for events with large |At| that are not 
modeled by either the signal or background PDFs, an outlier background Gaussian, Poh is 
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also included. The probability that an event is a signal event is assigned from the Mbc — AE 
PDF based on the event's position inside the signal box, 

/Sig IMbe, Ai., rj _ ^ ^^^^^^^ ^^^^^ ^ _ ^ ^^^^^^^^ ^^^^^ , [iZ) 

where p (r) is the purity in the signal region of each r-bin. Finally, the full time-dependent 
CP violation PDF is written as, 

P (At, q, Mbc, AE, r) = (1 - /oi)(/sig {M^c, AE, r)Psig (At, q) + 



i[l - /sig (Mbc, AE,r)]PBkg (At)| + i/oiPoi (At). 

(13) 

A control sample, D*^D^K~ , is selected for which no CP violation is expected. 

For this sample, the time- dependent fit reveals no CP violation, 

Acp (i?° ^ D*+D°K-) = 0.09 ± 0.13 (stat) . 
Sep (P° ^ D*+D^K-) = 0.01 ± 0.20 (stat), ^ ' 

where Sep represents the component of mixing induced CP violation. We then determine 
the CP asymmetry parameters of P° D*^ D*~ K^, 

Acp = -O.Olt°:l (stat) ± 0.09 (syst) 
P) sin 201 = OMtoil (stat) ± 0.06 (syst), ^ 

and the fit results are shown in Fig. [31 

The sources of systematic errors are addressed as follows. The uncertainties in the vertex 
reconstruction are the dominant sources of error. These include uncertainties in the IP 
profile, charged track selection based on track helix errors, helix parameter corrections. 
At and vertex goodness-of-fit selection, and SVD misalignment. Toy MC samples showed 
some small fitting bias for CP parameters due to low statistics in each sample. We take 
this bias as a systematic. The parameters for w and Aw, resolution function, physics 
parameters, background shape and signal fraction were varied by ±lcr. We estimate that 
the CP eigenstate component of the background does not exceed 10% of the total. The 
effect oi CP asymmetry in the background was estimated by changing the background PDF 
to be composed of 10% CP eigenstates with CP parameters at the physical limits. We also 
investigate the effect of mis-reconstructed signal events using signal MC. This is achieved 
by comparing the fit result of the signal MC sample with another fit on the same sample, 
which required that the events were reconstructed correctly. Tag-side interference comes 



from CP violation on the tag-side [17(], and is estimated with B D*lv. We generate MC 



pseudo-experiments and perform an ensemble test to obtain systematic biases [18|. The 
systematics are summarized in Table [TTl 

In the 3-parameter fit, the time-dependent At distribution assuming no direct CP viola- 
tion, is given by 

g-|At|/TgO c 

Psig(At,g) = -^——}^l-q/\w + q{l-2w)x 

r]y-^ cos {Anid At) — f^— ^ sin20i + r^^,— ^ cos20i^ sin (AmjAt) |. 
Jo \ Jq Jo J \ } 

(16) 
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FIG. 3: 2-parameter fit results of D*^ D*^ K^. (a) and (b) show the background subtracted 

At distribution for poor tags (0.0 < r < 0.5) and good tags (0.5 < r < 1.0). Background 
subtraction is performed by subtracting the background PDF from the raw At distribution. The 
solid curve represents tagged candidates and the dashed curve represents B^ tagged candidates, 
(c) and (d) show the B^ — B^ raw asymmetry, ( A^'^o — A^'^o ) / {N^o + N^o ) , where N^o {N^o ) is the 
number of B^ (B^) tags in At for poor tags and good tags, respectively. 



We obtain the time-dependent CP parameters, 

Jc/Jo = 0.601^:1 (stat) ± 0.08 (syst) 
2Jsi/Josin20i = -0.17ini(stat) ±0.09 (syst) 
2Js2/^ocos20i = -0.23l|]:g(stat) ±0.13 (syst). 



(17) 



and the corresponding fit curves are shown in Fig. HI 
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Resolution function 
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I—' rn 7"C7 1 /" * c T "\ o T* o vn mT" tHY* o 

jriiybicb paiaincLcio 


n nni 

U.UUl 


n nni 

U.UUi 


Background At 


0.016 


0.019 


Signal fraction 


0.005 


0.004 


Mis-reconstructed events 


0.014 


0.038 


Tag-side interference 


0.053 


0.001 


Total 


0.090 


0.058 



TABLE II: Time-dependent CP asymmetry systematics: 2-parameter fit. 

The systematics are determined in much the same way as in the 2-parameter fit. However, 
since BaBar claims evidence for a. ^ D^i (2536)D*^ contribution to the D*^ D*~ Kg 
final state [l^, and resonant structure is more important in the 3-parameter fit, we excluded 
the (2536) region, m {D*^K^) < 2.6 GeV/c^, and took the difference in the results as a 
systematic. The systematics are summarized in Table lllli 



Source 


S{Jc/Jo) 


(5(2J,i/Josin2(/>i) 


(5(2Js2/Jo cos 20i) 


Vertex reconstruction 


0.066 


0.071 


0.081 


Wrong tag fraction 


0.021 


0.033 


0.032 


Resolution function 


0.014 


0.039 


0.045 


Fit bias 


0.029 


0.018 


0.031 


Physics parameters 


0.005 


0.004 


0.003 


Background At 


0.005 


0.019 


0.012 


Signal fraction 


0.020 


0.009 


0.008 


Mis-reconstructed events 


0.020 


0.007 


0.026 


Tag-side interference 


0.000 


0.002 


0.004 


D+ (2536) 


0.019 


0.013 


0.069 


Total 


0.084 


0.094 


0.127 



TABLE III: Time-dependent CP asymmetry systematics: 3-parameter fit. 



In summary, we measure the branching fraction (Eq. E]) and CP parameters of 
D*+D*-K^ with a 2-parameter fit (Eq.[I5]), and 3-parameter fit (Eq.[I7D with 414 fb"^ at the 
Belle detector. An upper limit on the branching fraction of the intermediate two-body decay 
B^ — >■ D^i (2536)i5*^(Eq. [8]) is also obtained. In the 2-parameter fit, there is no evidence for 
direct CP violation and a large dilution of sin20i is implied. In the 3-parameter fit, Jc/Jo 
appears to be non-zero; a large value of this parameter could indicate the presence of a broad 
unknown D**~^ state. As in the 2-parameter fit, a large dilution of the CP asymmetry from 
polarization and resonant structure is implied since 2Jsi/ Jq sin 20i is consistent with zero. 
According to 2Js2/<^, is expected to be positive if this unknown wide resonance exists. 
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A t (ps) A t (ps) 




FIG. 4: 3-parameter fit results of D*^ D*^ Kg. (a) and (b) show the background subtracted 

At distribution for good tags (0.5 < r < 1.0) in the region, s~ < s"*" and > . The sohd curve 
represents tags and the dashed curve represents B^ tags, (c) and (d) show the B^ — B^ raw 
asymmetry in At for each half-Dahtz space and for good tags only. 



The parameter, 2Js2/<^ cos20i, has been measured, however, a model-dependent inference 
on the sign of cos 20i is not possible with current precision. 
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